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Abstract. In [5] the authors present an algorithm to parametrize approximately e- 
rational curves, and they show in 2 examples that the Hausdorff distance, w.r.t. to the 
Euclidean distance, between the input and output curves is small. In this paper, we ana- 
lyze this distance for a whole family of curves randomly generated and we automatize the 
strategy used in [5]. We find a reasonable upper bound of the Hausdorff distance between 
each input and output curve of the family. 



Introduction 

The problem of the approximate parametrization of algebraic plane curves goes as follows: 
given a plane algebraic curve C (that is the perturbation of a rational plane curve) and 
a tolerance e > 0, we want to find a new curve C, being rational, as well as a rational 
parametrization of it such that C and C are at certain small distance dependent on e. 

In [2] and [3] it was seen how to approximately parametrize algebraic plane curves and 
algebraic surfaces having an e-singularity of maximum multiplicity. In [5], using the tech- 
niques of e-points developed in [2] and [4], we extended the results in |2] to the general 
case of algebraic affine plane curves. More precisely, we provide in [5] the approximate 
parametrization algorithm which given C returns an approximate parametrization V(t) of 
the curve C. 

A natural question arises, the closeness between the input and output curves of the 
algorithm. In our case, this closeness notion is given by the Hausdorff distance. That is, we 
say that the input and output curves are close if their Hausdorff distance (as real curves) is 
small related to the tolerance. We proved in [5] that the Hausdorff distance between C and 
C is finite. 

In this paper, apply the approximate parametrization algorithm to a family of curves T 
of degree 4 which was randomly generated. For each curve C in T we compute a bound B 
of the Hausdorff distance between C and C. For all of the curves we get B < 2.2 and we 
obtain evidences that the actual distance is experimentally < 0.1. 

The paper is organized as follows. We recall the approximate parametrization algorithm 
in Section [TJ In Section [2] we explain how the family J- of curves of degree 4 was generated. 
The last section is devoted to the analysis of the distance between the the curves of T and 
their approximate parametrizations output by our algorithm. 

The following terminology will be used throughout the paper. || • || and || • ||2 denote the 
polynomial oo— norm and the usual unitary norm in C 2 , respectively. | • | denotes the module 
in C. The partial derivatives of a polynomial g € C[x, y] are denoted by := d t+ ^ g / d l xd^ y , 
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where v = 6 N 2 ; we assume that g° = g. Moreover, for it = £ N 2 , 117*1 = i + j. 
Also, it = (1, 0), and i| = (0, 1). 

1. Recalling the parametrization algorithm 

In this section we recall the algorithm presented in [5] as well as its main properties; see 
[5] for further details. We start with a fixed tolerance e, and with the implicit equation 
f(x,y) of a real plane algebraic curve C of exact degree d, which is the perturbation of a 
rational curve. C is supposed to satisfy that: 

(1) the degree d of / is proper; this means that there exists a partial derivative of /, of 
order d, that in module is strictly bigger than e||/||. 

(2) / is e-irreducible. 

(3) C has d different points at infinity, and it does not pass through (1 : : 0), (0 : 1 : 0). 

Conditions (1) and (2) guarantee that, under the tolerance, we really have an irreducible 
curve of degree d. Condition (3) ensures that C is either compact (as a subset of R 2 ) or it 
follows real asymptotes. Therefore, we are excluding curves having a parabolic behavior. 
The requirement on (1 : : 0), (0 : 1 : 0) is technical and it can be achieved by performing 
a suitable and orthogonal linear change of coordinates. 

The theoretical argumentation of the algorithm is as follows. First, the notion of exact 
singularity is replaced by the concept of e-singularity, similarly with the notions of exact 
multiplicity and e-multiplicity. Here, the first complication appears since the number of 
e-singularities is bigger than (expected) in the exact case; probably due to the perturbation. 
In order to deal with this difficulty we associate to each e-singularity a radius, and hence 
we see it as an Euclidean disk. Next, we introduce an equivalence relation on the set of 
disks and we define the e-singular clusters as the equivalence classes. Then, we define the 
e-multiplicity of the cluster as the maximum of the e- multiplicities within the class, and we 
take as canonical representant of the cluster an e-singularity where the e-multiplicity of the 
cluster is achieved. In this situation, we say that C is e-rational if the clusters satisfy the 
well-known genus formula of the exact case. More precisely, if {£luster ri (Qi)}i = i t ^ tS is the 
cluster decomposition (Qi denotes the canonical representant and rj the e-multiplicity of the 
cluster), we say that C is e-(affine) rational if 

s 

(d-l)(d-2)-^ri(ri-l) = 0. 

i=l 

Now, let us assume that C is e-rational, an let us see how the approximate parametrization 
algorithm proceeds. The basic idea is, as in the exact case (see [7]), to construct a suitable 
linear system of curves of degree d — 2. More precisely, if 

{€luster ri (Qi)}i=i t ..., s , Qi = (qn : q%2 ■ 1) 

is the cluster decomposition, we compute d — 3 simple e-points on C, say {Pi, . . . ,P^_3} 
with Pi = (pn : pi2 : 1). Again, we associate to each e-point an Euclidean disk via a radius, 
and we apply the equivalence relation. If, somehow, any Pj is identified with another Pj or 
with a singular cluster, we replace Pj by a new e-point. In this situation, we consider the 
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effective divisor 

s d—3 
i=l i=l 

and the (exact) linear system H of curves of degree (d — 2) given by D. That is, H is the 
linear system of curves of degree d — 2 having Qi as (rj — l)-base points, and Pi as simple base 
points. If we were working exactly, all intersection points in 7~L D C would be fixed (namely 
those points in D) with the exception of one point that would provide the parametrization. 
Indeed, in the exact case, the parametrization would be 

/ Res y (H(x,y, I), f) Res x (H(x, y, 1), f) 

vn-=i(* - a,!)'*'- 1 ) ntfc - ' uu(y - Qny^-v nti(y - 

being H(x,y,z) the homogeneous polynomial defining %. In the approximate case, instead 
of the exact division above, we take the quotient of the Euclidean division of each numerator 
by the corresponding denominator. 

The output curve derived from this process has the same structure at infinity as the input 
curve and the same degree, see [5], Theorem 4.5. These properties will play a fundamental 
role in the error analysis (see Section [3]). We outline the algorithm derived from the above 
ideas. 

(1) Compute the singular cluster decomposition {<£luster ri (Qi)}i=i,..., s ; say Qi = (q^i : 
%2 ■ !)• 

(2) If Y,U\ r i( r i - 1) 7^ - - 2), RETURN "C is not (affine) e-rational". If s = 1 
one may apply the algorithm in [2]. 

(3) Compute (d — 3) e-simple points {Pj}i<j<d-3 of C. Take the points over R, or as 
conjugate complex points. After each point computation check that it is not in the 
cluster of the others (including the clusters of Qi); if this fails take a new one. Say 
Pi = (Pi,l : Pi,2 ■ !)• 

(4) Determine the linear system T~L of degree (d — 2) given by the divisor Ylt=i r iQi + 
YltZfPi- Let H(t, x, y, z) = -ffi(x, y, z) + tH 2 (x, y, z) be its defining polynomial. 

(5) If [gcdCFfo^O), #!(*,!/,())) + 1] and [gcd(F(x, y, 0), H 2 (x, y, 0)) + 1] replace H 2 
by H 2 + p\x d ~ 2 + p 2 y d ~ 2 , where pi, p 2 are real and strictly smaller than e. Say that 
gcd(F(x, y, 0), H 2 (x, y, 0)) = 1; similarly in the other case. 

(6) Si(x,t) = Res y {H(x,y,l),f) and S 2 (y,t) = Res x (H(x, y, 1), /). 

(8) Ai(x) = YlUifr ~ %i) n{n - l) nS (* " 

My) = UU(y - ^Y^- l) n£f (» - 

(9) For i = 1, 2 compute the quotient l?j of by w.r.t. either x or y. 

(10) If the content of B\ w.r.t x or the content of B 2 w.r.t. y does depend on t, RETURN 
"degenerate case". 

(11) Determine the root Pi(t) of B±, as a polynomial in x, and the root p 2 (t) of i?2, as a 
polynomial in y. 

(12) RETURN V(t) = ^(t) ,p 2 (*)) . 
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2. Generating a family of 6-rational curves 



In this section, we generate the family of curves that will be used in the error analysis. 
We fix three points P x = (2 : : 1),P 2 = (0 : : 1) and P 3 = (1 : 1 : 1) in P 2 (C) and we 
consider the linear system of curves of degree 4 defined by the divisor 2Pi + 2P 2 + 2P3. Its 
defining polynomial is 

G(x,y,z,ui, ... ,u e ) = u 2 y 2 z 2 +u 3 y 3 z + u 4: y 4 + u 5 xyz 2 -(2u 2 + 3u 3 + Au 4l + ^u 5 + 2uQ)xy 2 z + 
u$xy 3 + uix 2 z 2 + (— tjU 5 + 2u 3 + 4u 4 + 2u 6 - u\)x 2 yz + (u 2 + u 3 + ^u 5 + \u\ + Ui)x 2 y 2 - 
U\x 3 z + {\u§ — u 3 — 2ii4 — U6 + \ui)x 3 y + \u\x A . 

Note that for every specialization of ui such that G(x,y, z,U\, . . . ,Uq) is irreducible, we 
get an (exact) rational curve. 

Now, for j = 1, . . . , 6 and i = 1, . . . , 10 let be a random integer number in the interval 
[0, 100]. We obtain 60 different polynomials Gij(x, y, z), j = 1, . . . , 6, i = 1, . . . , 10 setting 



in G(x,y,z,u\, . . . ,uq). Given i € {1,...,6} and j G {1, . . . , 10} we obtain a random 
perturbation gij(x,y) € ]R[x,y] of Gij(x,y,l) as follows 



where r\,r 2 , r 3 are integer numbers taken randomly in the interval [0, 100] and e = . The 
polynomials gij(x, y), j = 1, . . . , 6, i = 1, . . . , 10 have proper degree 4 and define 60 curves 
Cij verifying (1 : : 0), (0 : 1 : 0) ^ (C^j is the projective closure of Cy) and such that 
they have 4 different points at infinity. Therefore, each of 60 curves satisfies the hypothesis 
required in parametrization algorithm. 

Using the parametrization algorithm described in Section [H we conclude that 28 of the 
60 curves are e— rational. We show those curves in Fig. [TJ An statistical error analysis was 
given in 0. The precise equations of Cij as well as the parametrizations provided by the 
algorithm can be found in 

http:/ /www.aq.upm.es/Departamentos/Matematicas/srueda/fam4.pdf 



We describe the theoretical strategy for the error analysis that will be employed in the 

last section. Let e > be the tolerance, C the input curve and C the output curve provided 

by the approximate parametrization algorithm. Also, let f(x, y) and f(x, y) be the defining 

polynomials of C and C, respectively. Moreover, let V(t) be the parametrization of C output 

^ 

by the algorithm. Let C and C denote the real part of C and C, respectively. 



Now, we consider the Hausdorff distance H, induced by Euclidean distance d in the real 
plane. That is, for A, B C M 2 \ {0} 




k = j 
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9ij(x,y) = G ij {x,y,l) + €-^(x + y) + e 2 -^{x 2 + xy + y 2 )+ 
e 3 ^-(x 3 + x 2 y + xy 2 + y 3 ) 



3. Theoretical strategy for the error analysis 



H(A,B) = max{sup aeA {d(a, B)},sup beB {d(b, A)}} 
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FIGURE 1. Plot of the 28 e-rational curves Cjj randomly generated. 

where d(a,B) = inffe e ^{d(a, b)}, and by convection i?(0, 0) = and, for ^ A C X, 
H(A, 0) = 00 (see [T], for further details). 

The most relevant fact, for our analysis, is that H(C ,C ) < 00 (see [5], Lemma 6.1) 

TO IK 

because we want to bound H(C , C ). For this purpose, we will proceed in a similar manner 
to Section 6 in [5], not only for a couple of examples, but for the whole family of e-rational 
curves randomly generated in Section [2j More precisely, we consider the normal line to C 
at the generic point V(t), say Ci(t,s), as well as the normal line to C at the generic point 
(a, 6) G C K , say £2(0,, b,s). Moreover, we introduce the polynomials 

Di(M) = /(A(M))e5p)H, P 2 (a,M)=7(£ 2 (a,M)) eC(5)[*], 

where K(i) denotes the algebraic closure of R(i) and C(C) the field of rational functions over 
C. In addition, for every to G such that X , i(to,s) is well defined and has real roots, and 
for every (ao, 60) €= C R , such that T>2{ao, bo, s) is well defined and has real roots, we take 

pf{t ) = min{|s |/£>i(io,so) = and s G R}, 

pf (a , b ) = min{js | / T> 2 (a Q , b , s ) = and s G R}. 

Then, the supremum of ^(t), /of (a, b) provides an upper bound of the Hausdorff distance; 
at least for those subsets of both curves where the considered minimums are well defined. 

Because of computational difficulties, in our analysis, instead of computing p^(t), p\^(a, b), 
we will study 

pi(to) := min{|s | / £>i(*o, s o) = 0}, 
P2(ao,b ) := mm{\s \ / T> 2 (a , b , s ) = 0}. 

These quantities bound d(V(to),C) and d((ao,bo),C) respectively (here, we understand that 

d is the unitary distance in C 2 ), instead of d(V(to),C ) and d((ao,bo),C ). So, for those 
subsets of both curves, where the corresponding polynomials are well defined, we bound 

A(C R ,C R ) := max{su PteR {d(^(*),C)} > sup (0i6)eC E{d((a l 6),C)}}. 
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For this purpose, in the next section, we will follow the next steps: 

(1) We compute a bound of sup t£R {d(V (t) , C)} , by applying Corollary 6.2 in [5] that 
ensures that 

where Ai(t) ^ and 1 < i < n > , 





A (t) ' 


{(?) 


4(t) 



where D^t, s) = A n {t)s n H h A (t). 

(2) In order to bound sup ^u{d((a, b),C)}) one may apply the same corollary as 

in (1). However, it implies to maximize a bivariate rational function under the 
constrain f(x,y) = 0. This can be done, for instance, using Lagrange multipliers. 
Nevertheless, in practice, this is unfeasible. Instead, wc estimate the bound by taking 
a lattice of points (a, b) G C where we bound d((a, b),C). 

(3) Note that the quantity A(C , C ) gives information on how close every real point on 

_]g> 

each of the curves is of a complex point on the other curve. However, A(C , C ) < 
H(C ,C ). To avoid this difficulty, in our analysis, we will look for empirical ev- 
idences indicating that the computed bound of A(C , C ) also bounds H{C ,C ); 
for that we test empirically that, in our computations, pi (to) = P^(*o)i P2(oo, &o) = 
pf(a ,6 ). 

(4) It may happen for some to 6 M that (to) i s n °t well defined or simply that 
the bound it provides of d(V(to), C K ) is not satisfactory. Observe that to bound 
d(V(to), C K ) we can use the intersection of any line through V(to) with C. So in 
some cases we will also proceed in the following way. 

Let us consider the line at the generic point V(t) in the direction given by h € R 

n , \ (—, \ 2/i , , h 2 - 1 

Chit, S) = \ Pl[t) + ,P2\t) + 1 



We introduce the polynomial P/j(t,s) = f(Ch(t,s)) £ R(t,/i)[s]. For a fixed /io € 
R and for every to £ R, such that , D/j (to,s) is well defined and has real roots, 
d(V(t ),C R ) < pf (to), where 

p^(t ) = min{|s | /X>ho(to, so) = and s G R}. 

Thus, the supremum of (t) and pf (a, 6) provides an upper bound of the Hausdorff 
distance. Then for to G R such that D/ l() (to, s) is well defined, we may study 

Ph (to) :=mm{\so\/V ho (to,so) = 0} 

which is an upper bound of d(V(to),C) and plays the role of pi (to) i n the previous 
steps. 

4. Execution of the error analysis 

Let J- be the family of 28 e-rational curves of degree 4 defined in Section 3. For each curve 

C in J 7 we explain next how the bound of A(C , C ) was computed and show evidences that 



this bound is also an upper bound of the Hausdorff distance H(C ,C ). 
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4.1. Bound of sup t ^{d(V (t) , C)}. Let us denote by A the domain of T>i(tQ, s). If the curve 
is compact then the polynomial T>\(to, s) is well defined for every to £ R, otherwise it is not 
defined for two real poles /3±, P2 of V(t). In the family T only two curves are compact. 
In order to bound sup t£ ^{d(V (t) , C)} , we obtain an upper bound of pi(t) when t £ A. 

For this purpose we maximize the functions R±(t) = 4 

step (1) in Section [3]) in A as follows. 

Let a\ and 02 be the real roots of the denominator of R\{t). For all the curves in T , 
a\ and 02 are not real roots of the denominator of Ri{£). Let Ij, i = 1,2 be an interval 
isolating a, from aj, j 7^ i and from the real poles of Rz(t). Observe that Ri(t) and R%(t) 
are continuous in M \ (I\ U I 2) and in the adherence I\ U I2 of I\ U I2, respectively. We 
compute 

Bi = max{i?i(t) I i e R \ (Ji U h)} and B 2 = max{R 2 (t) \ t G hUl 2 }. 

Then B = max{i3i, #2} is an upper bound of p\{t) in A. 

The last column of the next table contains the computed bound B for each one of the 
curves Cj, i = 1, . . . , 28 of the family J 7 . 

Table 1: 



i 


01 


B 2 


B 


1 


0.3012751472 


1.784885546 


1.784885546 


2 


0.1680336313 


0.8228821157 


0.8228821157 


3 


0.2209183305 


1.143210796 


1.143210796 


4 


0.2457462218 


1.388890611 


1.388890611 


5 


0.4775061243 


1.471164469 


1.471164469 


6 


0.1854050321 


0.9172323537 


0.9172323537 


7 


0.3392516285 


1.238494405 


1.238494405 


8 


0.1687631697 


0.9278483955 


0.9278483955 


9 


0.4481254299 


1.345341665 


1.345341665 


10 


0.1706747632 


1.252669418 


1.252669418 


11 


0.4336254993 


1.328637472 


1.328637472 


12 


0.04502452088 


0.7028506083 


0.7028506083 


13 


0.2511290220 


1.849173820 


1.849173820 


14 


0.7973544750 


0.5426224779 


0.7973544750 


15 


1.947190823 


1.201605769 


1.947190823 


16 


0.1658993167 


2.124343900 


2.124343900 


17 


0.06346428265 


1.634020447 


1.634020447 


18 


1.401107905 


1.830395156 


1.830395156 


19 


1.690902532 


0.9731554792 


1.690902532 


20 


1.590285558 


1.659392056 


1.659392056 


21 


0.1543267485 


0.9761129297 


0.9761129297 


22 


0.8247063503 


0.8459442935 


0.8459442935 


23 


0.2452763324 


0.6140170288 


0.6140170288 


24 


0.08434729326 


0.7159251709 


0.7159251709 


25 


0.6464253153 


2.150679036 


2.150679036 



and R 2 (t) = (f) 



A (t) 
AW) 
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Table 1: continued 



26 


3.604620794 


1.418101314 


3.604620794 


27 


0.08779082555 


0.7809344831 


0.7809344831 


28 


0.3673519642 


1.773877016 


1.773877016 



We will improve next the bound given for curves C\% and C26- For a fixed ho € M we can 
write T>h (t, s) = B n (t)s n + • • • + Bo(t) to which Corollary 6.3 in j5] applies. Hence we can 
obtain an upper bound B h ° of sup tgK {(i('P(i), C)} maximizing the new functions 



B (t) 



Bi(t) 



and (2) 



B (t) 



B 2 (t) 



in the domain of T>h (t,s). As described earlier for Ri(t) and R2(t) we obtain respectively 



and B, 

For curves Ciq and C26 we computed the upper bound B ho of sup tm {d(V (t) , C)} using 
different values of ho and we found bounds improving the ones given earlier for the values 
of ho shown in the next table. 

Table 2: 



ho 



% 


ho 


Bf° 




B h 


16 


-1 


0.5870746534 


1.287063889 


1.287063889 


26 


1 

20 


0.2525792337 


0.001139009266 


1.178706930 



4.2. Empirical bound of sup, b ^^u{d((a, b), C)}). In this section we estimate the bound 

of sup -,m{d((a, b),C)}). We estimate the bound by taking a lattice of points (a, b) £ C 

where we bound d((a,b),C) estimating p^(a,b). We show evidences for p^(a,b) being small 
and for p2{a, b) = p^(a, b). 

If the curve is not compact, first we analyze the behavior of the input and output curves 
through the real asymptotes. Let us suppose that C is a non compact curve in T and let 
C\ and £2 be its real asymptotes. By [5], Corollary 4.6 the real asymptotes of C and C are 
parallel lines so the Hausdorff distance between them can be easily computed. Let C\ and 
£2 be the real asymptotes of C parallel to L\ and £2 respectively. We the value of 

7 7 = max{iJ(£i,£i) ) iJ(£ 2 ,£ 2 ))} 

for all the non compact curves of T in the next table. Then we proceed as follows: 

(1) For each negative integer i we compute the set fij of intersections of C R with the line 
x = i. We obtain := max{p^( a > b) \ (a, b) € Qi} and m; := max{p2(o 5 b) \ (a, b) E 
Oj}, and we check that = mj. 

(2) We repeat the previous step until 

mm{\ P 1{i, b) - F(£i,£i)|, \pf (i, b) - #(£ 2 ,£ 2 )| / (t, b) e n t } < e. 

(3) Let T\ be the smallest value of i until termination of this process. 

We perform this experiment also for each positive integer i to obtain in this case the 
highest value t% such that the inequality in step (2) holds. At the same time we check that 
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mf~ = rrii with positive integers i = 1, . . . ,t 2 - The same process is repeated for y = j, to 
obtain the negative and positive integers T3,T4, respectively such that 

minilpf (a,j) - H(d,£i)\, \/$(a,j) - H(C 2 ,C 2 )\/ (aj) € Q?} < e 

where Qi is the set of intersections of C K with the line y = j. Let (m- J ) R := max-fp^ (a, b) | (a, b) 6 

f^} and m? : = max{p 2 (a, 6) | (a, 6) € f2 J }. We also check that (m- 7 ) = m? with integers j = 
r 3 , . . . , -1, 1, . . . ,r 4 . Let ?n := max{m i ,m J ' [ i = n, . . . , -1, 1, . . . ,r 2 , j = r 3 , . . . , -1, 1, . . . ,r 4 }. 

Let [ri,T2] x [t3,T4]. We empirically consider that out of the compact [ri,T2] X [t3,T4], 
the curves behave as the asymptotes, and the empirical bound of sup, -^u{d((a, b),C)}) 

in [tl,t 2 ] x [r 3 ,r 4 ] is m. 

The next table shows the compact set [ti,T2] x [t3,T4] obtained for e = 10~ 6 in all the 
curves except for numbers 25 and 28 for which we took e = 1CT 5 . The reason being that 
we run out of memory before reaching the box outside of which the curves behaved like the 
asymptotes with e = 10 -5 . 

Table 3: 



l 


T l,T 2 


U3, 


r 4. 


m 


V 


1 


[-3434, 3428] 


[-2069, 


2066] 


0.04474051996 


0.002685992105 


2 


[-6730, 6732J 


[-3833, 


3835] 


0.01909150476 


0.007250422655 


3 


[-2120, 2120] 


[-1030, 


1031] 


0.02523781400 


0.0008098244306 


4 


[-2485, 2492] 


[-1370, 


1374] 


0.03061351675 


0.0006391265474 


5 


[-4157, 4160] 


[-2000, 


2003] 


0.03559861599 


0.0007910063013 


6 


[-11115, 11089] 


[-4663, 


4642] 


0.02312094507 


0.005224665954 


7 


[-22154, 22121] 


[-9370, 


9334] 


0.03942395261 


0.01044227346 


8 


[-11207, 11214] 


[-4756, 


4761] 


0.02118572337 


0.005399010161 


9 


[-28777, 28684] 


[-11887, 


11825] 


0.05222449790 


0.01308589690 


10 


[-1922, 1925] 


[-1063, 


1065] 


0.02633886946 


0.0003845822760 


11 


[-4555, 4552] 


[-1912, 


1912] 


0.03226327617 


0.001373344655 


12 


[-5017, 4997] 


[-2114, 


2107] 


0.01246202377 


0.002415851115 


14 


[-797, 780] 


[-242, 


235] 


0.01907109331 


0.001277331144 


15 


[-4809, 4774] 


[-1410, 


1396] 


0.07529126612 


0.005987329671 


16 


[-3841, 3866] 


[-15344, 


15340] 


0.05528214429 


0.01084165622 


17 


[-139, 16] 


[-955, 


832] 


0.03702449872 


0.008212462120 


19 


[-13558, 13560] 


[-6428, 


6446] 


0.03068560995 


0.01035036823 


20 


[-4290, 4291] 


[-2287, 


2288] 


0.04051580912 


0.001538264208 


21 


[-14465, 14520] 


[-5275, 


5300] 


0.01835471004 


0.005065273865 


22 


[-2286, 2282] 


[-982, 


982] 


0.01911087028 


0.0002824836230 


23 


[-946, 940] 


[-400, 


399] 


0.01154559037 


0.0004413541720 


24 


[-157, 3724] 


[-1215, 


1213] 


0.01530235590 


0.0001629862393 


25 


[-18910, 18839] 


[-5882, 


5863] 


0.08880952924 


0.03464418857 


26 


[-3400, 3398] 


[1769, 


1769] 


0.03324362713 


0.001139009266 


27 


[-308, 5274] 


[1659, 


1642] 


0.01648328102 


0.0002351747177 


28 


[-6279, 6210] 


[-2001, 


1982] 


0.03876376237 


0.01250853150 
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If the curve C is compact we consider a compact set [ti,T2] X [73, T4] containing C . Then 
we compute m as previously described checking also that mf = mi,(m^) = m? with 
i = n, -1, 1, ■ ■ ■ ,r 2 , j = r 3j ...,-l, 1, ... ,r 4 . 



Table 4: 



i 


[n,r 2 ] 


N,r 4 ] 


m 


13 
18 


[-9/512,8041/1024] 
[-5/512,592745/512] 


[-20057/1024,2117/1024] 
[-8723/1024, 304847/1024] 


0.04595703645 
0.09228397972 



4.3. Empirical evidences. Now, we perform some empirical tests to show evidences that 
pf(t) is smaller than the upper bound B of sup tgK {d(V (t), C)} given in Section [4.1i First, 
let -Di(s) = limt,_ 5 .-|- 00 Pi(t, s). Then, for every curve of the family J- let 

X = min{|so| / -Di(so) = and so G R}. 

We checked that x = mm {|so| / ^i( s o) — 0} in all cases. Since the roots of a polynomial 
depend continuously on its coefficients, for every 5 > there exists K > such that for 
all |to| > K there is a root sq of T>i(to,s) with \\x — S0II2 < It may happen that these 
roots are all complex. However, in our example, we see that pi(( — W) k ) = Pi({— W) k ) for 
k = 1, . . . , 20. Let xi = mm{pf ((-10) fe ) \k = 1, . . . , 20} and X 2 = max{pf ((-10) fc ) | k = 
1, . . . , 20}. We show these computations in the next table. 



Table 5: 



i 


X 


Xi 


X2 


1 


0.001918863706 


0.001918568088 


0.001922644324 


2 


0.004169957700 


0.004161970065 


0.004170583666 


3 


0.0006994105148 


0.0006993543405 


0.0007001275847 


4 


0.0006662568567 


0.0006659169716 


0.0006706114665 


5 


0.001185494963 


0.001184862374 


0.001191850748 


6 


0.0002278441391 


0.0002266218453 


0.0002434582060 


7 


0.004927221514 


0.004922425041 


0.004975309877 


8 


0.001581018062 


0.001580991493 


0.001581020714 


9 


0.003893937744 


0.003759724534 


0.008421643120 


10 


0.0004038643651 


0.0004036487366 


0.0004066244124 


11 


0.001324473283 


0.001324364325 


0.001325864354 


12 


0.0006513574356 


0.0006500882228 


0.0006514569249 


13 


0.004018853976 


0.004018487265 


0.004023539023 


14 


0.001399552125 


0.001399552100 


0.001399552395 


15 


0.006540790926 


0.006540789722 


0.006540806320 


16 


0.01873734902 


0.001829670103 


0.03745550824 


17 


0.0003018645266 


0.0002989553843 


0.0003020916294 


18 


0.003558391244 


0.003557440039 


0.003570581899 
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Table 5: continued 



i 


X 


Xi 


X2 


19 


0.0005206627618 


0.0005145010748 


0.0005996180488 


20 


0.001782196531 


0.001779734571 


0.001782429636 


21 


0.001157562963 


0.001155811228 


0.001179946145 


22 


0.001096140685 


0.001095297129 


0.001096220492 


23 


0.0004670320359 


0.0004670115956 


0.0004670339681 


24 


0.001057066460 


0.001056738383 


0.001061264741 


25 


0.06749573507 


0.005148188392 


0.06814019159 


26 


0.001564435599 


0.001563775730 


0.001572903801 


27 


0.001268468809 


0.001267989563 


0.001274605947 


28 


0.0005550274941 


0.0004502065315 


0.0005632086626 



Observe that for most of the curves of the family T the bound B computed in Section 5.1 
was reached at B 2 . Let us suppose that the isolating intervals of the real poles of R\(t) are 
h = [in, ^12] and I 2 = ^21^22]- Then if T 2 = {t € R / R' 2 {t) = 0} then 

B 2 = max{R 2 (t)/t G T 2 U {t n , t 12 , t 21 , t 22 }}. 

Furthermore, only for curve C25 the set T 2 is nonempty. We show in the next table 
that v = m&x{pf (tn), pf{t\ 2 ), pf (t 2 \), pf(t 22 )} is very small compared to the value of 
/i = m&x{R 2 (tii), R 2 (t\ 2 ), R 2 (t 2 i), R 2 (t 22 )} and also we compare them with B 2 . 

Table 6: 



i 


B 2 


/" 


V 


1 


1.784885546 


1.784885546 


0.01895037003 


2 


0.8228821157 


0.8228821157 


0.003187256654 


3 


1.143210796 


1.143210796 


0.007932988669 


4 


1.388890611 


1.388890611 


0.01141905470 


5 


1.471164469 


1.471164469 


0.01307915412 


6 


0.9172323537 


0.9172323537 


0.004427051081 


7 


1.238494405 


1.238494405 


0.008019188716 


8 


0.9278483955 


0.9278483955 


0.004590148599 


9 


1.345341665 


1.345341665 


0.009661551521 


10 


1.252669418 


1.252669418 


0.009307873805 


11 


1.328637472 


1.328637472 


0.01086359027 


12 


0.7028506083 


0.7028506083 


0.002603616572 


13 


1.849173820 


1.849173820 


0.01407639891 


14 


0.5426224779 


0.5426224779 


0.002434019507 


15 


1.201605769 


1.201605769 


0.01147127202 


16 


2.124343900 


2.124343900 


0.01663838395 


17 


1.634020447 


1.634020447 


0.01149183389 


18 


1.830395156 


1.830395156 


0.01586113294 


19 


0.9731554792 


0.9731554792 


0.006071538946 


20 


1.659392056 


1.659392056 


0.01644172109 


21 


0.9761129297 


0.9761129297 


0.005852977427 


22 


0.8459442935 


0.8459442935 


0.004712176306 


23 


0.6140170288 


0.6140170288 


0.002448022081 
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Table 6: continued 



i 


B 2 




V 


24 


0.7159251709 


0.7159251709 


0.006089518478 


25 


2.150679036 


2.150679036 


0.05229161121 


26 


1.418101314 


1.418101314 


0.01218396938 


27 


0.7809344831 


0.7809344831 


0.007553005576 


28 


1.773877016 


1.773877016 


0.01971360528 



For curves C i: i = 14, 15, 19, 26 the bound B is achieved in B\. Let T 1 = {t G R / R[(t) = 
0}, then 

Si = max{{fli(t)/i € Ti U {in, t 12 , t 2 i, ^22}} U {lim t ^ ±OOJ Ri(t)}}. 

In fact, in those curves Si is equal to max{i?i(i) / t G Ti}. 

For each real pole of Ri(t) as well as for each real critical value of R\(t), we consider a 
sequence of isolating intervals of length l/10 fc+5 , we take the middle point tk, and we 
analyze pi(ifc), (ifc). After a certain k\ the sequences become stable, let Efa be the set 
containing /^(tfa) for the feith element of each one of the sequences constructed for each 
real pole of R\{t). Let 71 = maxi?^. Similarly, after a certain k 2 the sequences for the real 
critical values of R\(t) become stable. Let be the set containing pf(ifc 2 ) for the foth 
element of each one of the sequences constructed for each real critical values of Ri{£). We 
call 72 = maxEk 2 - If the curve is not compact, we also perform this experiment for each of 
the two real poles /3i,/?2 of V(t). The sequences to /3i,/?2 become stable after a certain 
and we get a set E^ 3 and 73 = max E^ ■ 



Table 7: 



i 


7i 


72 


73 


1 


0.05760637790 


0.07103885930 


0.002685991017 


2 


0.03680694646 


0.04030133537 


0.007250422354 


3 


0.03985346560 


0.05219935882 


0.0008098246652 


4 


0.04358984318 


0.06256427563 


0.0006391277723 


5 


0.05222536651 


max{0.009579373594, s Cl } 


0.0007910068620 


6 


0.03997031881 


0.03468777768 


0.005224667630 


7 


0.05420521510 


0.06230995426 


0.01044227515 


8 


0.03842889832 


0.03215152707 


0.005399009187 


9 


0.06425437603 


0.08158624352 


0.01308589929 


10 


0.03857067352 


0.03800183984 


0.0003845830022 


11 


0.05050891041 


max{0. 007217929868, s C5 } 


0.001373344987 


12 


0.02590356741 


0.01124398891 


0.002415851766 


13 


0.07620545140 


0.05039734679 




14 


0.02750848667 


max{0.01846235849, s C5 } 


0.001277332295 


15 


0.06016762755 


0.2971598318 (*) 


0.005987330302 


16 


0.09714515451 


0.03765565867 


0.01084165864 
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Table 7: continued 



17 





06901149454 


0.01524196447 




0.008212460400 


18 





07011950234 


max{0.01367185130, s C5 


} 




19 





05269958334 


0.2771889722 (*) 




0.01035036672 


20 





05744093868 


max{0.01461432871,s C2 


} 


0.001538264277 


21 





05357323512 


0.03486848807 




0.005065270849 


22 





04568079074 


max{0.009770701627, s C2 } 


0.0002824834209 


23 





02825492360 


0.06513903382 




0.0004413527932 


24 





05459990275 


0.01907914217 




0.0001629875241 


25 





09628252969 


0.1341804080 




0.03464419053 


26 





05149574516 


1.431046152 (*) 




0.001139009002 


27 





06133149851 


0.01990495615 




0.0002351734340 


28 





09919600166 


0.09135201774 




0.01250853076 



4.4. Behaviour of pf(t). In the grate majority of our computations Pi(t) is denned and 
Pi(t) = pi(t). In some cases for a given to S R the set {|so| / fi(to> s o) = and sq € M} is 
empty or pf(t) and pi(t) happen to be different. Then it should be taken into consideration 
that there exists ho £ R such that p^ (to) < oo. 

Let s c denote {p^(tfe)} for the sequence {tk} to the critical point c of Ri(t). For curves 
number 5, 11, 18, 20 and 22 the sequence s c could not be computed for some critical point 
c of R\(t). For example, curve C5 has 6 critical points and the sequence s ci for critical 
point ci was not defined. We write max{0.009579373594, s Cl } where 0.009579373594 is the 
maximum of the values at which the rest of the sequences stabilized. 

The results for 72 marked with (*) indicate that Px~(t k ) 7^ Pi(t) f° r the sequence {t k } to 
one of the critical points of Ri(t). For curve number 26 there are 6 critical values of R\(t) 
and 1.357539211 < \pf(t k ) - p\(t k )\ < 1.357539223 for the sequence {t k } to the critical 
value C5. 

For each one of the curves highlighted and for the critical point of R\(t) where p^(t) did 
not behave properly (there was only one of those points of each curve) we proceed as follows. 
We consider a sequence of isolating intervals J k of length l/10 fc+5 , we take the middle point 
t k , and we analyze p^ (t k ) for different values of ho- The next table shows the values of ho 
giving good results together with the value p^it^) at which the sequence became stable 
and the new maximum -y' 2 - 



Table 8: 



i 


ho 


J* (+ \ 


7 2 


5 


1 


0.05275311956 


0.05275311956 


11 


1 


0.05120027918 


0.05120027918 


14 


1 


0.06065146651 


0.06065146651 


15 


1 


0.09855121223 


0.09855121223 


18 


1 


0.08937284288 


0.08937284288 


19 


4/5 


0.1099404739 


0.1099404739 


20 


1 


0.07070538339 


0.07070538339 


22 


1 


0.06827542251 


0.06827542251 


26 


1 


0.07787356026 


0.07787356026 
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